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Abstract 

The Kinetic Gas theory iike two-agent money exchange model, recently introduced 
in the Econophysics of wealth distributions, is revisited. The emergence of Boltzmann- 
Gibbs like distribution of money to Pareto's law in the tail of the distribution is exam- 
ined in terms of 2 x 2 transition matrix with a general and simplified outlook. Some 
additional interesting results are also reported. 
Introduction 

Econophysics of Wealth distributions [1] is an emerging area where some Statistical 
Physicists and Economists have been engaged in interpreting real economic data of 
money, wealth and income of all kinds of people pertaining to different societies and 
nations. Economic activities have been assumed to be analogous to elastic scattering 
processes [3, 4, 5]. Analogy is drawn between Money (m) and Energy (E) where 
temperature (T) is average money (< m >) of any individual at equilibrium. Some 
early attempts [2] have been made to understand the income distributions which follow 
Pareto's law (P(m) oc l/m a+1 ) at the tail of the distributions with the index (a) 
varying from 1 to 2.5. 

Kinetic Gas theory like models are introduced to exploit apparent similarities be- 
tween a many particle system and a social system of many agents. The primary attempt 
is to understand the distributions of money/ income corresponding to different classes 
of people. We revisit such a Kinetic theory like two-agent money exchange model as 
recently proposed by Chakrabarti and group [4, 6]. In this model any two agents chosen 
randomly from a number of agents (A) are allowed to interact (trade) stochastically 
and thus money is exchanged. Stochasticity is introduced in terms of a parameter 
< e < 1 into the interaction. One arrives at a Boltzmann-Gibbs (exponential)-type 
distribution (P(m) oc exp(— (5m)) of individual money. In the next stage, a saving 
propensity factor (A) is incorporated to show that the distribution shifts away from 
the exponential distribution. A peak appears at a value other than at zero. In the 
later stage the saving propensity factor is made random (among the agents) but frozen 
in time. This brings a distribution qualitatively different from the earlier: one gets 
a power law at the tail of the distribution. This indicates the emergence of Pareto's 
Law (P(m) oc l/m a+1 ) at the tail of the distribution which is well known among 
Economists/ Econophysicists. In the above class of models total money (M = J2i m i) 
of all the agents is invariant in time. Also the money (m) is conserved locally which 
means the sum of money of two agents before and after trade (interaction) remains 
constant: mi(t + 1) + mj(t + 1) = mi(t) + mj(t). 
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We systematically examine the role played by the parameters e and A in this con- 
served class of models and attempt to understand the emergence of behaviours in 
terms of probability distribution functions {P{m)) of money (m). Also we come across 
additional interesting features in the model with variable saving propensity (A). 

The time evolution of money in the model without saving propensity is described 
as follows: 

m i (t+l) = e(m i (t)+m j (t)) (1) 



mj (t + 1) = (1 - e)(mi(t) + mj(t)), (2) 

Therefore, we can say that the distribution in money is evolved through 2x2 
transition matrices (T): 

r =(lle ll e ) 

The above matrix is, however, singular which means the inverse of this matrix does 
not exit. This indicates that the evolution through such transtion matrices is bound 
to be irreversible. As a result we get exponential (Boltzmann) distribution of money. 
This can be perceived here in a different way too. When we take a product of such 
matrices and we get back one which is nothing but the leftmost matrix: 





The above signifies the fact that when two agents happen to interact repeatedly 
(via this kind of transition matrices), the last of the interactions is what matters (the 
last matrix of the product survives). This 'loss of memory' may be attributed to the 
path to irreversibility in time here. 

Let us play the game in a different way. Suppose two agents trade (interact) in 
such a manner that we may arrive at the following general transition matrix: 

Tl= (i-61 1-62 ) 

where e\ and 62 are two independent random fractions varying between and 1 (but 
not equal). This simply signifies that e\ fraction of money of the 1st agent added with 
€2 fraction of money of the 2nd agent is retained by the 1st agent after the trade. The 
rest of their total money is shared by the 2nd agent. This is a little generalization of the 
earlier case. However, this matrix is now nonsingular (as long as e\ 7^ €2) and the two- 
agent interaction process may be said to be reversible in time. Therefore, we expect to 
have an equilibrium distribution of money which may be qualitatively different from 
the earlier exponential (Boltzmann-Gibbs like) one. In the above general matrix, if we 

put €i = l and €2 = 0, this reduces to an Identity matrix / = ( ^ ) which is the 

stationary and trivial case (no interaction). 

In Fig.l we plot two distributions (for e\ = €2 and e\ 7= £2) in support of the above 
discussions. The computer simulation is performed on a system of iV=1000 agents. 
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Two agents are selected randomly to interact (trade). No qualitative change is seen 
in the distributions when we take different system sizes. A single interaction between 
any two agents is defined here as a single time step. Simulation is done for 10 5 time 
steps and averaging is done over 10 4 initial configurations. The distributions, we plot 
here, are not normalized. 




Figure 1: Distribution of money (m) in two cases: Exponential curve is for the case when the 
transition matrix is singular (no time reversal symmetry), the other curve with a peak is for 
the case when the transition matrix is made non-singular by choosing e x ^ e 2 . 

Let us now look at the case where a saving propensity factor A is incorporated [4]. 
When two traders/ agents meet, each of them is supposed to keep aside A- fraction 
(fixed) of their individual money (This is like a gamble with precaution.). Rest of 
the money ((1 — A)-fraction of the sum of their total money) is redistributed (between 
them) with the stochasticity factor e (0 < e < 1). The transition matrix now looks 
like: 

( A + e(l-A) e(l-A) \ 

^ (l-e)(l-A) A + (l-e)(l-A) ) 

We may now rescale the matrix elements by assuming e\ = A + e(l — A) and 
62 = e(l — A) in the above matrix. Therefore, the above transition matrix reduces to 

T2= (i-61 1-62 J 

Thus the matrix T2 is of the same form as T\ . Here, as < e < 1 and A is something 
between and 1, we do have < ei < 1 and < £2 < 1. As long as e\ and £2 are 
dfferent, the determinant (A = e\ — €2 = A) of the matrix is nonzero. Therefore, 
the effect of the saving propensity factor A, thus introduced, essentially leads to have 
non-singular transition matrix. Hence it is clear from the above discussion that the 
distribution (in money) would likely to be qualitatively no different from what can 
be achieved with transition matrices like T\ with e\ and €2 (ei 7^ €2), in general. 
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The distributions obtained for different A (as reported in [6]) may correspond to the 
difference in ei and €2 as we see. 

In the next stage, when the saving propensity factor A is made random, the tran- 
sition matrix between any two agents having different A's (say, Ai and A2) would now 
look like: 



Again we rescale the elements by putting ei = Ai + e(l — Ai) and £2 = e(l — A2). Hence 
this matrix can also be reduced to the same form as that of T2. 



The determinant here is A = e± — £2 = Ai(l — e) + eA2- Here also A is ensured to be 
nonzero as all the parameters e, Ai and A2 have the same range: between and 1. This 
means that each transition matrix for two-agent money exchange remains non-singular 
which signifies the interaction process to be reversible in time. We may also check here 
that < ?i < 1 and < £2 < 1, which again corresponds to the transition matrix T\ as 
discussed before. Therefore, it may be apparent that qualitatively different distributions 
are possible when we appropriately tune the two elements e\ and €2 in the general form 
of transition matrix T\ (We have not done so here in this paper to explicitly demonstrate 
that further.). Nevertheless the incorporation of the parameter A goes a step closer to 
interpret real economic data. However, the emergence of power law tail (Pareto's law) 
in the distribution is not well understood in the model[6] we are discussing. 

In this context, we present an additional interesting feature of the model (as pro- 
posed in [6]) with variable saving propensity (A). Suppose, we consider A to have only 
two fixed values Ai and A2 and that they are widely different. This may be thought 
of the society to have only two kinds of people: some of them do save a very large 
fraction (fixed) of their money and the other kind of people who save a very small 
fraction (fixed) of their money. Introducing this binary-A in the model brings out an 
interesting feature: a double-peak distribution in money which can be seen from Fig. 2 
(This result is also obtained by simulating a system of ^=1000 agents.). Thus two dis- 
tinct economic classes appear out of this. It has been observed that the system evolves 
towards a distinct two-peak distribution as the difference in Ai and A2 is increased 
systematically. Later it is seen that we still get two-peak distributions even when Ai 
and A2 (one high and the other low) are distributed in a narrow range around their 
fixed values. 

The class of models here we discuss is conserved. The sum of money of all the 
agents M = J2i m i is fixed for all time. This is ensured by the two-body interaction 
process (rule of the game) where the sum of money of the two agents is conserved before 
and after interaction as mentiond in the beginning. So the sum of two elements of a 

column of a transition matrix T = ( | n | 12 ) has to be unity by design: t\\ + £21 = 

1, £12 + £22 = 1- Whatever extra parameter we add in the model, no matter, the 
matrix has to retain this property. (However, an extra parameter, like A, may help 




Ai + e(l-Ai) e(l-A 2 ) 
(l-e)(l-Ai) A 2 + (l-e)(l-A 2 ) 
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Figure 2: Double-peak distribution of money(m) with fixed values, Ai=0.2 and A2=0.8: 
Emergence of two Economic classes. 

explaining things better.) Therefore, this kind of conserved models, in general, may be 
understood in terms of a general transition matrix like T\ as discussed in the beginning. 
T\ essentially tells that one agent retains randomly ei-fraction (0 < e\ < 1) of his own 
money added with random e2-fraction (0 < 62 < 1) of the other, where the other keeps 
the rest. 
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